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This  paper  presents  a  comparative  study  of  twelve  equivalent  circuit  models  for  Li-ion  batteries.  These 
twelve  models  were  selected  from  state-of-the-art  lumped  models  reported  in  the  literature.  The  test 
data  used  is  obtained  from  a  battery  test  system  with  a  climate  chamber.  The  test  schedule  is  designed 
to  measure  key  cell  attributes  under  highly  dynamical  excitations.  The  datasets  were  collected  from  two 
types  of  Li-ion  cells  under  three  different  temperatures.  The  multi-swarm  particle  swarm  optimization 
algorithm  is  used  to  identify  the  optimal  model  parameters  for  the  two  types  of  Li-ion  cells.  The  usefulness 
of  these  models  is  then  studied  through  a  comprehensive  evaluation  by  examining  model  complexity, 
model  accuracy,  and  robustness  of  the  model  by  applying  the  model  to  datasets  obtained  from  other  cells 
of  the  same  chemistry  type. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Hybrid  electric  vehicles  (HEVs),  plug-in  hybrid  electric  vehi¬ 
cles  (PHEVs)  and  battery  electric  vehicles  (BEVs)  are  being  actively 
developed  and  deployed  to  achieve  significant  fuel  consumption 
and  carbon  emission  reductions  in  many  markets  throughout  the 
world.  Traction  battery  packs  are  presently  the  performance  and 
cost  bottlenecks  of  these  electrified  vehicles.  To  ensure  safe,  reli¬ 
able  and  efficient  operations  of  the  traction  batteries  under  the 
most  demanding  and  grueling  driving  conditions,  an  effective  bat¬ 
tery  management  system  (BMS)  must  be  used.  A  key  function  of 
the  BMS  is  to  monitor  the  conditions  and  states  of  the  traction 
battery  pack,  such  as  State  of  Charge  (SOC)  and  State-of-Health 
(SOH).  Since  these  variables  are  not  directly  measurable  by  any 
sensors,  they  need  to  be  inferred,  commonly  from  model-based 
estimation  algorithms  [1-7].  Therefore,  accurate  battery  models 
are  of  utmost  importance.  In  addition  to  providing  accurate  estima¬ 
tions,  it  is  important  to  strike  a  balance  between  model  complexity 
and  accuracy  so  that  the  models  can  be  embedded  in  microproces¬ 
sors  and  provide  accurate  results  in  real-time  [8].  In  other  words, 
it  is  important  to  have  models  that  are  accurate  enough,  and  not 
unnecessarily  complicated. 

Electrochemical  models  that  aim  to  capture  all  key  behaviors  of 
battery  cells  [9-1 1  ]  often  can  achieve  high  accuracy.  They  are  suit¬ 
able  for  understanding  the  distributed  electrochemistry  reactions 
in  the  electrodes  and  electrolyte  [8].  However,  they  typically  deploy 
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partial  differential  equations  with  a  large  number  of  unknown 
parameters.  The  complexity  often  leads  to  significant  requirement 
for  memory  and  computation.  In  addition,  they  frequently  run  into 
over-fitting  problems  due  to  their  significant  number  of  param¬ 
eters  and  the  model  robustness  under  extrapolation  is  typically 
poor— rendering  their  practical  usefulness  questionable.  Hence, 
PDE-based  electrochemical  models  are  not  desirable  for  actual  bat¬ 
tery  management  in  electrified  vehicles— if  lumped  battery  models 
are  accurate  enough,  they  are  likely  to  be  the  preferred  choice. 

Equivalent  circuit  battery  models  have  been  studied  especially 
for  the  purpose  of  vehicle  power  management  control  and  bat¬ 
tery  management  system  development.  They  are  lumped  models 
with  relatively  few  number  of  parameters.  For  example,  Plett 
[2]  proposed  several  lumped  models,  including  the  simple,  zero- 
hysteresis,  one-state  hysteresis,  combined,  and  enhanced  self 
correcting  (ESC)  models.  In  [2],  test  datasets  of  an  HEV  Li-polymer 
cell  was  used  to  compare  these  models.  The  resistance-capacitance 
(RC)  network  based  equivalent  circuit  models  were  also  widely 
studied,  such  as  the  first-order  RC  [12-14],  second-order  RC 
[15-17]  and  third-order  RC  models  [18].  Dedicated  elements 
depicting  the  battery  hysteresis  behavior  were  sometimes  added 
to  the  RC  models,  such  as  the  first-order  RC  model  with  hystere¬ 
sis  [19-21]  and  third-order  RC  model  with  hysteresis  [22],  etc.  In 
the  literature,  there  were  few  studies  comparing  these  commonly 
used  equivalent  battery  models.  The  goal  of  this  paper  is  to  system¬ 
atically  compare  the  practicality  of  these  models,  including  model 
complexity,  model  accuracy  under  both  training  and  validation 
data,  and  generalizability  of  these  models  to  multiple  cells.  The 
last  attribute  is  important  because  with  improvement  in  battery 
chemistry,  manufacturing  process,  cooling  systems,  and  cell-to-cell 
balancing  circuitry,  it  is  likely  the  cell  variations  will  continue  to 
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Arbin  BT  2000  Tester 


Fig.  1.  Configuration  of  the  battery  test  bench. 


decrease,  to  such  an  extent  that  we  may  feel  comfortable  (or  be 
driven  economically)  to  measure  and  monitor  the  temperature  and 
voltage  of  every  few  cells  instead  of  each  and  every  cell.  To  real¬ 
ize  this  vision,  it  is  important  to  confirm  that  a  model  calibrated 
based  on  the  measurements  of  one  cell  can  be  applied  to  predict 
the  behavior  of  several  adjacent  cells  placed  under  the  same  load 
and  are  operating  under  similar  conditions  such  as  environmental 
temperature. 

In  this  paper,  we  study  twelve  commonly  used  lumped  bat¬ 
tery  models.  They  are  compared  using  multiple  cell  datasets  we 
acquired  under  different  temperatures  for  two  types  of  Li-ion 
cells.  The  multi-swarm  particle  swarm  optimization  (MPSO)  [23] 
is  implemented  to  identify  the  optimal  model  parameters  based  on 
the  training  data  from  a  single  cell  (arbitrarily  selected  from  eight 
cells).  Then,  the  usefulness  of  these  models  is  compared  by  using 
validation  test  data  obtained  from  eight  cells  of  the  same  chemistry. 
This  process  was  repeated  for  two  types  of  Li-ion  batteries. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section 
2,  the  experimental  setup  to  acquire  the  data  is  described.  In  Sec¬ 
tion  3,  the  twelve  equivalent  circuit  battery  models  are  described. 
The  MPSO-based  model  identification  process  is  introduced  in  Sec¬ 
tion  4.  The  model  comparison  results  are  discussed  in  Section  5, 
followed  by  conclusions  presented  in  Section  6. 


2.  Experimental  setup 

2A.  Battery  testing  systems 

The  experimental  setup  is  shown  in  Fig.  1.  It  includes  an 
Arbin  BT2000  tester,  a  thermal  chamber  for  environment  con¬ 
trol,  a  computer  for  user-machine  interface  and  data  storage,  a 
switch  board  for  cable  connection,  and  the  battery  cells.  During  the 
charging/discharging,  voltage,  current,  temperature  of  each  cell  is 
measured  and  recorded  at  10  Hz. 

An  impedance  measurement  system  (see  Fig.  2)  was  also 
designed  to  record  the  battery  impedance  behavior  at  various  exci¬ 
tation  frequencies.  It  comprises  a  data  acquisition  unit  for  signal 
generation  and  data  collection,  a  power  amplifier  for  signal  ampli¬ 
fication,  and  two  shunt  resistors  for  current  measurement.  The 
measurement  of  impedance  is  obtained  from  sinusoidal  excita¬ 
tions,  in  which  both  voltage  and  current  in  sinusoidal  forms  (with 
dc  offset)  are  recorded  and  their  complex  quotient  is  computed  as 
the  cell  impedance.  To  reduce  measurement  noise,  a  convolution- 
based  method  is  adopted  to  identify  the  phase  delay  of  impedance 
using  the  Parseval’s  Theorem  for  magnitude  calculation  [24].  Only 
two  cells  from  each  Li-ion  cell  types  are  tested  and  the  sinusoidal 
input  is  modulated  on  top  of  a  precisely  generated  dc  offset  to 
ensure  we  are  not  charging  or  discharging  the  cells  during  the 


Table  1 

Main  specifications  of  the  LiNMC  and  LiFeP04  cells. 


Type 

Nominal 

Nominal 

Upper  cut-off 

Lower  cut-off 

capacity  (Ah) 

voltage  (V) 

voltage  (V) 

voltage  (V) 

LiNMC 

0.94 

3.70 

4.20 

2.50 

LiFeP04 

1.10 

3.30 

3.60 

2.00 

test.  The  sampling  time  is  input  frequency  dependent  and  can  vary 
between  1  / 40  and  5  kHz. 

2.2.  Battery  test  schedule 

Two  types  of  cylindrical  Li-ion  cells  are  selected  for  our  tests. 
One  is  lithium  nickel-manganese-cobalt  oxide  (LiNMC)  UR14650P 
cells  from  Sanyo  and  the  other  is  lithium  iron  phosphate  (LiFeP04) 
APR18650M1A  cells  from  A123.  Eight  cells  each  were  purchased 
on  the  open  market.  Their  key  specifications  are  shown  in  Table  1. 
These  cells  were  placed  in  cell  holders  in  the  thermal  chamber.  They 
are  independently  tested  using  1 6  channels  of  the  battery  tester.  For 
the  cells  of  the  same  chemistry,  the  same  loading  profile  is  applied. 

The  test  schedules  shown  in  Fig.  3  are  designed  to  generate  rich 
excitations  for  the  two  types  of  cells.  Each  experimental  procedure 
takes  about  two  weeks  to  finish.  In  each  procedure,  it  begins  with 
a  characterization  test  at  temperature  T=10°C,  followed  by  two 
identical  tests  conducted  at  T=  35  °C  and  T=  22  °C.  A  static  capac¬ 
ity  test,  a  hybrid  pulse  test,  a  dc  resistance  test,  a  DST  test  and 
a  FUDS  test  are  consecutively  conducted  in  each  characterization 
test.  The  purpose  of  the  static  capacity  test  is  to  measure  the  cell 
capacity.  The  hybrid  pulse  test  is  a  sequence  of  pulse  cycles.  Each 
pulse  cycle  is  composed  of  the  standard  Hybrid  Pulse  Power  Charac¬ 
terization  (HPPC)  profile  and  a  self-designed  discharging/charging 
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Fig.  3.  Flow  chart  of  the  test  schedule. 
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pulse  profile.  The  self-designed  profile  is  a  combination  of  pulses 
with  different  amplitudes  and  durations.  It  serves  two  functions: 
to  move  the  cell  SOC,  and  to  excite  the  cell  dynamically.  The  self- 
designed  pulse  profile  for  LiNMC  cells  is  shown  in  Fig.  4.  That  for 
LiFeP04  cells  is  similar.  The  pulse  amplitudes  are  merely  enlarged  to 
achieve  the  same  current  rates  as  those  for  LiNMC  cells.  The  dc  resis¬ 
tance  test  uses  the  standard  testing  program  from  Arbin  to  estimate 
the  internal  resistance.  The  Dynamic  Stress  Test  (DST)  and  Federal 
Urban  Dynamic  Schedule  (FUDS)  tests  excite  the  cells  in  driving 
cycle-based  conditions.  Between  two  adjacent  test  points  in  each 
characterization  test  the  cells  are  charged  or  discharged  to  reach 
the  desirable  initial  SOC  values  (around  90%)  and  rested  to  reach 
cell  equilibrium.  After  these  tests,  the  impedance  test  is  conducted 
to  collect  data  for  battery  impedance  spectrum  analysis.  Then,  the 
aging  cycles  are  conducted  at  T=  22  °C.  In  each  aging  cycle,  the  cell 
is  charged  or  discharged  at  a  constant  rate  until  the  cut-off  voltages. 

2.3.  Datasets  for  this  work 

The  datasets  collected  in  the  three  characterization  tests  for 
fresh  cells  in  the  first  experimental  procedure  (before  aging  cycles) 
are  used  for  the  model  identification  and  comparison  in  this  paper. 
The  hybrid  pulse  test  datasets  under  three  different  temperatures 
for  the  reference  cells  (LiNMC  cell  in  Channel  17,  LiFeP04  cell 
in  Channel  25)  are  used  as  the  training  datasets.  The  DST  and 
FUDS  datasets  for  the  reference  cells  serve  as  the  model  validation 
datasets.  In  addition,  the  hybrid  pulse  test,  DST  and  FUDS  datasets 
under  three  different  temperatures  for  all  the  sixteen  cells  are  used 
to  assess  the  generalization  of  the  optimized  models  to  multiple 
cells.  Since  the  SOC  of  the  battery  packs  in  electrified  vehicles  is  not 
allowed  to  be  less  than  a  certain  threshold  (e.g.,  20%)  due  to  bat¬ 
tery  life  consideration,  the  portion  of  the  test  data  below  10%  SOC 
in  these  datasets  are  not  used  in  model  comparison.  The  current, 
voltage  and  SOC  of  the  reference  LiNMC  cell  in  the  training  dataset 
at  T=  22  °C  are  shown  in  Fig.  5.  Note  the  dataset  below  1 0%  SOC  was 
discarded.  Similar  pre-processing  was  done  for  the  DST  and  FUDS 
datasets  (see  Figs.  6  and  7,  a  portion  is  showed  to  give  more  details). 

3.  Model  structures 

A  total  of  twelve  equivalent  circuit  battery  models  are  selected 
from  the  literature  and  used  in  the  study.  These  models  were 
selected  to  form  a  comprehensive  subset  aiming  to  cover  the 
majority  of  the  lumped  model  templates  studied  in  the  past. 
The  effects  of  battery  SOC  and  current  are  covered  explicitly  in 
these  models.  In  addition,  the  hysteresis  observed  during  dynamic 
loading  (mainly  due  to  the  Li-ion  uneven  distribution  limited  by 
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Fig.  5.  Current,  voltage  and  SOC  of  the  reference  LiNMC  cell  in  the  training  dataset 
at  22  °C. 
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Fig.  6.  Current,  voltage  and  SOC  of  the  reference  LiNMC  cell  in  a  portion  of  the  DST 
dataset  at22°C. 
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Fig.  7.  Current,  voltage  and  SOC  of  the  reference  LiNMC  cell  in  a  portion  of  the  FUDS 
dataset  at22°C. 
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Table  2 

Twelve  battery  models. 

Model 

(1)  The  combined  model  [2] 


(2)  The  simple  model  [2] 


(3)  The  zero-state  hysteresis 
model  [2] 


(4)  The  one-state  hysteresis 
model  [2] 


(5)  The  Enhanced 
Self-correcting  (ESC)  model 
(two-state  low-pass  filter) 
[2] 


(6)  The  Enhanced 
Self-correcting  (ESC)  model 
(four-state  low-pass  filter) 
[2] 


(7)  The  first-order  RC  model 
[12-14] 

(8)  The  first-order  RC  model 
with  one-state  hysteresis 
[19-21] 

(9)  The  second-order  RC  model 
[15-17] 

(10)  The  second-order  RC 
model  with  one-state 
hysteresis 

(11)  The  third-order  RC  model 
[18] 

(12)  The  third-order  RC  model 
with  one-state  hysteresis 
[22] 


Model  equations  and  description 


Vk  =  /Co  -  Kt /zk  -  I<2zk  + 1<3  In (zk)  +  IQ  ln(l  -  zk)  -  R0Ik 

where  Vk,  zk  and  Ik  are  the  battery  terminal  voltage,  SOC  and  current,  respectively;  K0,  IQ,  IQ,  IQ  and  IQ  are  parameters 
for  the  dependence  between  open  circuit  voltage  (OCV)  and  SOC;  Ro  is  the  internal  ohmic  resistance  which  depends 
on  the  current  direction.  The  optimization  variable  vector  6  =  [IQ,  IQ ,  IQ,  IQ,  IQ,  R+,R~]. 

Vk  =  OCV(Zk)  -  RoIk 

where  OCV(Zk)  represents  the  dependence  between  OCV  and  SOC  in  the  form  of  a  table.  0  =  [Voc  i , . . . ,  Voc  u,  Rj,  Rq  ]• 

f  '*>*’ 

Vk  =  OCV(zk)-RoIk -skM,  sk  =  <  -1,  Ik  < -s, 
l^k-l,  \IK\  <  S, 

where  M  is  assumed  to  be  a  constant  coefficient  depicting  the  hysteresis  level  and  s  is  adequately  small  and  positive. 

6=  [Voc,  t,...,  Voc,  u,R+,R0,M]. 

hk+ 1  =  exp(-|/c4At|)hfc  +  [1  -  exp(-|/c4At|)]H 


Vk  =  OCV(Zfc)  -  R04  +  hk 

where  hk  is  the  hysteresis  voltage,  k  is  a  decaying  factor,  At  is  the  sampling  time  and  H  is  the  maximum  amount  of 
hysteresis  voltage  which  is  positive  for  charge  and  negative  for  discharge.  6  =  [Voc,i , . . . ,  Voc,i2,  Rq,R~,k,H+,H~ ]. 
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where /i  and  f2  are  two  states  of  a  low-pass  filter  of  the  current;  ai  and  a2  are  the  diagonal  elements  of  the 
state-transition  matrix  of  the  filter,  gi  and  g2  are  the  elements  of  the  output  matrix  of  the  filter,  and  a  zero-dc  gain 
constraint  is  applied  to  them.  6  =  [V0c,i  ,  y0c,i2 ,R%,R0,K,H+,H-,ai,a2,g\]. 
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Vk  =  OCV(Zfc)  -  Ro4  +  hk  +gl/l,k  +g2f2,k  +g3f3,k  +  g4f4,k 

where  a  four-state  filter  is  used.  0  =  [Voc,i, . . . ,  V0C)i2,  Rj,  R~,  k,  H+,  H~ ,ol\,ol2,  0L3,  014,  g\,  g2,  g3]. 
Ui,k+i  =exp(-At/ri)t/l!k+R1[l  -exp(-At/ri)]4 


Vk  =  OCV(zk)-R0Ik-U,,k 

where  U 1  and  X\  =  RiCi  are  the  voltage  and  time  constant  of  the  RC  network.  6  =  [Voc  ,i,  •  •  • ,  Voc,  12, R£,R0,R\,  ti]. 
Vk  =  OCV(Zk)  -  R04  -  U\,k+hk 

where  0  =  [VOC, 1 , . . . ,  Voc,12 ,  R+ ,  RQ  ,  k,  H+ ,  H~ ,  Ri ,  n  ] . 


Vk  =  OCV(zk)-R0Ik-Uu-U2,k 

where  U2  and  r2  =R2C2  are  the  voltage  and  time  constant  of  the  second  RC  network. 
0  —  [V0C,1,  •  •  •  ,  Voc, 12,  ^0  ’  ^0  ’  Tl’  ^2,  T2]. 

Vk  =  OCV(Zk)  -  Ro4  -  uu  -  U2M  +  hk 

where  0  =  [VocA, ...,  V0C,12,  R+,R~,k,  H+,  H~,Ru  n  ,R2,  r2]. 


Vfc  =  OCV(Zk)  -  RoIk  -  Uu  -  U2M  -  U3,k 

where  U3  and  t3  =R2C3  are  the  voltage  and  time  constant  of  the  third  RC  network. 

0  —  [  Voc,  1  ?  •  •  •  ,  Voc, 12,  ^0  ’  ^0  ’  Tl’  ^2,  X2 ,  R3,  T3]. 

Vk  =  OCV(Zk)  -  RoIk  -  uu  -  U2M  -  U3,k  +  hk 

where  0=  [Voc,  1, ...,  Voc,  12,  R+ ,  R0 ,  k,  H+ ,  H~ ,  Ru  R2,  t2,  R3,  r3]. 


diffusion)  is  included  explicitly.  The  effect  of  temperature  how¬ 
ever  is  not  captured.  The  effect  of  which  can  be  captured  through 
temperature-dependent  model  parameters.  These  battery  models 
are  summarized  in  Table  2. 


4.  Identification  of  optimal  model  parameters 

4. 1.  Optimization  variables  for  the  battery  models 

Out  of  the  twelve  models,  except  for  model  1,  the  OCV  is  mod¬ 
eled  as  a  function  of  SOC,  and  V0C,jJ  =  1,2,...,  12  is  used  to 
represent  the  jth  optimization  variable  for  describing  OCV.  When 
simulating  the  battery  model,  OCV  is  calculated  by  look-up  table. 
The  model  variables  to  be  optimized  for  all  the  twelve  battery  mod¬ 
els  are  all  summarized  in  Table  2.  In  this  table,  R  J  and  R~  are  internal 
ohmic  resistances  for  the  discharge  and  charge,  respectively;  H+ 


and  H~  denote  the  maximum  hysteresis  voltages  for  the  discharge 
and  charge,  respectively.  The  bounds  of  the  parameter  vector  to 
identify,  Oi  and  Gy,  were  specified  based  on  prior  knowledge  of  the 
two  types  of  Li-ion  batteries.  The  bounds  of  the  same  parameters  for 
different  models  are  kept  the  same  for  fair  comparison.  For  exam¬ 
ple,  the  upper  bounds  of  R+  in  all  the  models  were  set  to  0.2  £2  and 
0.1  £2  for  LiNMC  and  LiFePOq  cells,  respectively.  It  is  important  to 
note  that  \af  <  1,  j=  1,  . . .,  4  is  imposed  to  ensure  stability  of  the 
low-pass  filters  for  the  ESC  models. 

4.2.  Optimization  algorithm 

To  reduce  the  probability  of  trapped  at  a  local  minimum,  the 
PSO-based  global  optimization  approach  was  adopted.  Compared 
to  genetic  algorithm,  another  global  optimization  method  used  for 
battery  modeling  [16,22],  PSO  has  fewer  parameters,  has  a  more 
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Table  3 

The  numerical  process  of  the  multi-swarm  particle  swarm  optimization  (MPSO)  method  for  battery  model  optimization. 

Step  1 :  The  bounds  of  the  parameter  vector,  0L  and  0Ut  are  assigned  and  the  generation  index  gen  is  set  to  0. 

Step  2:  Within  the  space  bound,  randomly  generate  an  initial  swarm  (population)  P0  which  comprises  N  particles,  i.e.,  . . . ,  .  Each  particle  represents  a 

candidate  solution  of  the  model  parameters. 

Step  3:  Calculate  F(P0),  and  record  personal  bests  of  particles,  i.e.,  pbest°, . . . ,  pbest ° ,  where  F(  )  is  the  objective  function,  half  the  sum  of  the  squared  voltage  errors 
of  the  chosen  model  for  the  training  dataset. 

Step  4:  If  gen  is  less  than  the  maximum  generation,  Mgen,  perform  the  following  steps  (1 )— (3). 

(1)  Split  the  swarm.  Firstly  sort  the  particles  in  a  descending  order  according  to  their  superiority  levels  determined  by  the  objective  function  values.  Then,  the  first 
particle  is  chosen  as  the  local  best  of  the  first  sub-swarm.  Ns  -  1  particles  with  the  largest  Euclidean  distance  from  the  local  best  are  assigned  to  be  the  other 
members  of  the  first  sub-swarm.  Repeat  this  selection  procedure  for  the  remaining  particles  until  the  remaining  particles  are  less  than  Ns  so  as  to  establish  all  the 
sub-swarms. 

(2)  Evolve  each  particle  in  each  sub-swarm  by  taking  advantage  of  the  velocity  and  position  update  equations  of  Krohling  and  Coelho’s  PSO.  Note  that  the  local  best 
in  each  sub-swarm  is  updated  with  a  probability  of  85%  to  enhance  the  global  search  capability.  The  updated  particle  that  violates  the  boundary  constraint  is 
adjusted  to  meet  the  boundary  constraint. 

(3)  gen  =gen  + 1.  Calculate  F(  )  for  the  updated  swarm,  and  record  updated  personal  best  of  particles  according  to  their  superiority  levels,  i.e.,  pbest8™, . . . ,  pbest8™. 
Step  5:  If  the  generation  index  gen  is  equal  to  the  maximum  generation,  Mgen,  end  the  swarm  evolution.  Record  the  best  particle  in  the  final  swarm  and  the  best 

objective  function  value.  The  best  particle  is  the  optimized  model  parameter  vector. 


effective  memory  utilization,  and  is  more  efficient  in  maintaining 
the  diversity  of  the  search  candidates  [25].  A  hybrid  multi-swarm 
particle  swarm  optimization  (HMPSO)  algorithm  for  constrained 
optimization  was  developed  to  improve  the  global  search  capabil¬ 
ity  of  the  standard  PSO  in  [23  ].  In  the  new  algorithm,  the  differential 
evolution  (DE)-based  mutation  and  swarm  splitting  concept  are 
incorporated  into  the  standard  PSO.  Since  the  DE-based  mutation  is 
too  time-consuming  for  battery  modeling  on  large  datasets,  MPSO 
was  applied  without  DE.  For  MPSO,  at  each  generation,  the  swarm 
is  split  into  several  sub-swarms  and  each  sub-swarm  evolves  inde¬ 
pendently  by  taking  advantage  of  the  search  equations  of  Krohling 
and  Coelho’s  PSO  [26],  resulting  in  very  effective  diversity  of  the 
swarm.  The  general  framework  of  MPSO  for  the  battery  model 
parameterization  is  shown  in  Table  3.  The  readers  are  referred  to 
[23]  for  more  details  of  the  MPSO  algorithm. 

There  are  only  three  parameters  to  adjust  for  the  MPSO  algo¬ 
rithm,  shown  in  Table  4.  These  values  were  found  to  generate 
good  optimization  results  in  our  cases.  A  large  value  was  assigned 
to  the  maximum  generation  variable,  Mgen  when  identifying 
all  the  twelve  models  at  the  expense  of  slow  computations. 
For  each  model,  the  optimization  algorithm  was  individually 
executed  on  the  training  datasets  under  three  different  temper¬ 
atures  so  that  the  temperature  dependence  can  be  added  to  the 
model. 

5.  Model  optimization  results 


RC  models  improves  the  average  RMS  error  by  about  7.9%  com¬ 
pared  to  the  first-order  RC  model  without  hysteresis.  Adding  RC 
networks  can  also  improve  the  model  accuracy  however  adding 
more  complexity  beyond  2RC  network  is  not  helpful. 


5.1.  Model  comparison  using  the  training  and  validation  datasets 

The  objective  function  for  the  model  parameter  optimization 
was  defined  for  model  accuracy,  measured  by  the  average  root- 
mean-squared  (RMS)  error  between  the  test  datasets  and  the 
output  from  the  optimized  models.  Experimental  results  from  all 
three  temperatures  are  used.  However,  the  model  parameters  are 
allowed  to  change  and  be  optimized  for  the  three  temperatures. 

5.1.1.  LiNMC  battery  reference  cell 

The  maximum,  minimum  and  average  RMS  errors  of  the  twelve 
models  optimized  for  the  training  datasets  (hybrid  pulse  tests 
under  three  different  temperatures)  of  the  reference  LiNMC  cell 
are  shown  in  Fig.  8.  The  model  number  is  as  shown  in  Table  2.  It 
can  be  seen  that  the  RC  models  are  consistently  more  accurate  than 
the  first  six  models.  In  addition,  adding  one-state  hysteresis  to  the 


Table  4 

The  three  key  parameters  of  the  MPSO  algorithm. 

Parameter  Mgen  N  Ns 

Value  10,000  60  8 
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Fig.  9.  Model  validation  performance  for  the  reference  LiNMC  cell  under  the  DST 
datasets. 


364 


X.  Hu  et  al.  /  Journal  of  Power  Sources  198  (2012)  359-367 


25 


| 

8  20|- 

o 
> 

Lj 

5  is* 


<0 

i  io 
=> 

UL 


LU  5 
CO 
2 
cr 


- 1 - 1 - 1 - 1 - 1 - 

0 

— 1 — 

- 1 - 1 - 1 - 1 - 

O  Maximum  RMSE 

O  Average  RMSE 

*  Minimum  RMSE 

>  0 

Q 

o 

U 

O 

O 

- 

■  0  "  =  c 

The  first-order  RC  model 

* 

ft\ 

O 

O 

o 

o 

*  *  *  * 

0 

°  a  a  □  ( 

1 

_ 1 _ 1 _ 1 _ 1 _ 1 _ 

w 

_ 1 _ 

♦  .  .  *  . 

_ i _ i _ i _ i _ 

2  3  4  5  6  7  8  9  10  11  12 

Model  structure  number 


1 


o'* 

CL 


LU 

CO 

cr 


Model  structure  number 


Fig.  10.  Model  validation  performance  for  the  reference  LiNMC  cell  under  the  FUDS  F*S- 12.  Model  training  performance  for  the  reference  LiFeP04  cell, 

datasets. 


It  is  important  to  compare  the  model  performance  under  unseen 
validation  datasets,  using  which  the  model  prediction  capabili¬ 
ties  are  better  evaluated.  The  RMS  errors  of  the  twelve  models 
optimized  for  the  hybrid  pulse  tests,  but  evaluated  using  the  DST 
datasets  (validation  datasets)  under  three  different  temperatures 
are  shown  in  Fig.  9.  It  can  be  seen  that  all  twelve  models  have  worse 
performance  than  under  the  training  datasets.  Among  all  twelve 
models,  the  zero-state  hysteresis  model  has  the  worst  validation 
result,  indicating  that  its  robustness  might  be  a  concern.  Similar 
to  the  training  results  shown  in  Fig.  8,  the  RC  models  show  con¬ 
sistently  better  performance  compared  to  the  first  six  models.  The 
first-order  RC  model  is  almost  as  good  as  more  complex  models 
under  the  validation  datasets-which  seems  to  suggest  it  is  a  good 
choice  balancing  between  model  robustness  and  complexity. 

The  RMS  errors  of  the  twelve  models  under  the  FUDS  datasets 
(another  validation  datasets)  under  three  different  temperatures 


(a)  One  hybrid  pulse  cycle 
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Fig.  11.  Voltage  responses  of  the  first-order  RC  model  and  the  reference  LiNMC  cell 
under  22  °C,  (a)  one  hybrid  pulse  cycle  (training);  (b)  one  DST  cycle  (validation)  and 
(c)  one  FUDS  cycle  (validation). 


are  shown  in  Fig.  1 0.  Since  the  FUDS  test  is  a  driving  cycle-based  test, 
the  results  are  important.  Again  the  first-order  RC  model  stands  out 
as  the  best  model,  achieving  good  prediction  accuracy  while  using 
a  simple  model  structure. 

From  the  validations  results  above,  we  see  that  while  more  com¬ 
plicated  models  can  achieve  better  accuracy  under  the  training 
datasets,  they  seem  to  be  more  complicated  than  necessary,  and 
their  over-fitting  characteristics  does  not  help  when  the  model  is 
exposed  to  validation  datasets  that  it  has  not  seen  before.  Consid¬ 
ering  both  model  accuracy  and  model  robustness,  the  first-order 
RC  model  is  determined  to  be  the  best  model  for  the  LiNMC  cell. 
The  voltage  responses  of  the  first-order  RC  model  and  the  refer¬ 
ence  LiNMC  cell  for  one  cycle  in  each  dataset  are  shown  in  Fig.  11. 
The  optimized  battery  model  predicts  cell  voltage  accurately. 

While  the  numerical  results  clearly  show  the  superior  perfor¬ 
mance  of  the  first-order  RC  model,  one  cannot  but  ask  the  following 
question:  the  low-pass  filter  in  the  ESC  model  seems  to  be  con¬ 
ceptually  similar  to  the  function  of  a  single  RC  network.  So  why 
is  not  the  more  complex  model  achieving  better  results  under  the 
training  dataset?  The  first  reason  is  numerical  problem.  Despite 


Fig.  13.  Model  validation  performance  for  the  reference  LiFeP04  cell  under  the  DST 
datasets. 
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Fig.  14.  Model  validation  performance  for  the  reference  LiFeP04  cell  under  the  FUDS 
datasets. 


> 


Fig.  16.  Model  generalization  performance  under  the  hybrid  pulse  test  datasets  for 
eight  LiNMC  cells. 


the  fact  we  are  using  an  optimization  method  with  globally  sit¬ 
uated  particle  as  search  candidates,  the  solution  is  more  prone 
to  be  trapped  into  local  minima  for  more  complicated  models. 
Take  the  two-state  low-pass  filter  in  the  ESC  model  as  an  exam¬ 
ple.  Constraints  \aj\ <  1,  j=  1,  2(stability)  are  too  loose.  From  the 
optimization  result  of  the  second-order  RC  model  with  one-state 
hysteresis  and  a\  oc  e— At/^i  ci ,  a2  oc  e~At/R2C 2  ?  We  know  that  ,  a2 
are  in  the  intervals  of  [0.95,  1]  and  [0.98,  1].  (It  is  extremely  dif¬ 
ficult  to  know  the  two  tight  boundary  constraints.  Typically,  we 
just  know  \oij\ <  1,  j=  1,  2  from  our  prior  knowledge.)  After  tighten¬ 
ing  the  constraints  on  ot\  and  a2 ,  the  optimization  results  can  be 
improved.  From  numerical  optimization  viewpoint,  the  RC  models 
are  therefore  better  than  the  ESC  models  in  terms  of  insensitiveness 
to  boundary  constraints  and  avoiding  local  minima.  Although  the 
optimization  results  can  be  improved  by  tightening  the  boundary 
constraints  on  and  a2 ,  it  is  still  worse  than  that  of  the  second- 
order  RC  model  with  one-state  hysteresis.  This  may  be  due  to  the 
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Fig.  15.  The  voltage  responses  of  the  first-order  RC  model  with  one-state  hysteresis 
and  the  reference  LiFeP04  cell  under  22  °C,  (a)  one  hybrid  pulse  cycle  (training);  (b) 
one  DST  cycle  (validation)  and  (c)  one  FUDS  cycle  (validation). 


zero-dc  gain  constraint  for  the  ESC  models  which  does  not  apply  to 
the  RC  networks. 

5.1.2.  LiFeP04  battery  reference  cell 

The  RMS  errors  of  the  twelve  models  optimized  for  the  training 
datasets  (the  hybrid  pulse  power  tests  under  three  different  tem¬ 
peratures)  of  the  reference  LiFeP04  cell  are  shown  in  Fig.  12.  The 
RC-network  based  models  again  show  better  performance  com¬ 
pared  to  the  first  six  models.  It  also  seems  that  adding  a  one-state 
hysteresis  to  the  RC-network  models  is  helpful. 

Again  the  optimized  models  are  assessed  using  the  DST  and 
FUDS  cycles  to  validate  the  robustness  of  the  obtained  models.  The 
RMS  errors  of  the  twelve  models  using  the  DST  and  FUDS  datasets 
under  three  different  temperatures  are  shown  in  Figs.  13  and  14.  It 
can  be  seen  that  the  first-order  RC  model  with  one-state  hysteresis 
is  preferred  for  the  LiFeP04  cell,  due  to  its  balance  between  model 
complexity  and  accuracy.  The  voltage  responses  of  the  first-order 
RC  model  with  one-state  hysteresis  and  the  reference  cell  for  one 
cycle  in  each  subtest  in  the  characterization  test  under  22  °C  are 
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Fig.  17.  Model  generalization  performance  under  the  DST  datasets  for  eight  LiNMC 
cells. 
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cells. 


Fig.  20.  Model  generalization  performance  under  the  DST  datasets  for  eight  LiFeP04 
cells. 


shown  in  Fig.  15.  Again  the  model  prediction  output  is  accurate 
even  under  validation  datasets. 

5.2.  Generalization  to  datasets  from  multiple  cells 

Traction  battery  pack  in  an  electrified  vehicle  often  consists  of 
hundreds  or  even  thousands  of  cells.  Since  safety  is  a  major  concern 
for  Li-ion  batteries,  cell-level  management  and  control  are  very 
important.  For  battery  management,  we  can  monitor  and  model 
each  cell  with  its  own  sensors  and  models  at  the  expense  of  sig¬ 
nificant  hardware  and  computation  cost.  This  arrangement  is  used 
in  some  early  electrified  vehicles  because  of  the  large  cell-to-cell 
variations,  caused  by  manufacturing  variability,  thermal  gradient 
and  other  factors.  With  improved  battery  chemistry,  manufactur¬ 
ing  process  controls,  improved  cooling  systems,  and  cell-to-cell 
balancing  circuitry,  it  is  likely  the  cell  variations  will  continue 
to  decrease,  to  such  an  extent  that  we  will  feel  comfortable  to 
measure  and  monitor  the  temperature  and  voltage  of  every  few 


cells  instead  of  each  and  every  cell.  Under  this  envisioned  future 
scenario,  modeling  each  cell  is  unnecessary.  Instead,  it  is  possible 
we  will  generalize  the  model  established  and  adapted  for  a  sin¬ 
gle  cell  and  use  it  to  predict  the  behavior  of  several  adjacent  cells. 
Therefore,  it  is  important  to  see  which  model  structure  works  well 
when  the  model  trained  for  one  cell  works  on  measured  datasets 
from  other  cells.  Here,  the  models  obtained  from  the  hybrid  pulse 
testing  datasets  of  the  reference  cells  (LiNMC  cell  in  Channel  17, 
LiFeP04  cell  in  Channel  25)  were  used  to  predict  the  behavior  of 
other  cells  of  the  same  chemistry.  The  ability  to  “extrapolate  to 
multiple  cells”  is  measured  by  the  maximum,  minimum,  and  aver¬ 
age  RMS  errors  using  the  datasets  from  all  eight  cells  of  the  same 
electrochemistry  under  the  three  different  temperatures.  The  stan¬ 
dard  deviation  of  the  RMS  error  is  also  calculated  to  demonstrate 
the  spread  of  model  error  among  eight  cells. 

5.2.1.  LiNMC  cell 

The  RMS  errors  of  the  twelve  models  under  the  hybrid  pulse  test 
datasets  of  all  eight  LiNMC  cells  are  shown  in  Fig.  1 6,  all  of  which  are 
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Fig.  19.  Model  generalization  performance  under  the  hybrid  pulse  test  datasets  for 
eight  LiFeP04  cells. 


Fig.  21.  Model  generalization  performance  under  the  FUDS  datasets  for  eight 
LiFeP04  cells. 
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a  little  higher  than  those  in  Fig.  8,  as  one  would  expect.  The  aver¬ 
age,  min  and  max  RMS  errors  again  show  similar  trend.  Adding  RC 
networks  or  one-state  hysteresis  to  the  first-order  RC  model 
improves  the  model  generalization  capability  a  little.  The  max¬ 
imum  improvement  on  the  average  RMS  error  is  approximately 
2.8%.  The  RMS  errors  of  the  twelve  models  under  the  DST  datasets 
of  all  eight  LiNMC  cells  are  shown  in  Fig.  17.  Those  under 
the  FUDS  datasets  are  shown  in  Fig.  18.  It  can  be  seen  that 
the  first-order  RC  model  achieves  good  accuracy  with  minimum 
complexity. 

5.2.2.  LiFeP04  cell 

Similar  to  the  LiNMC  cells,  results  for  the  hybrid  pulse  test, 
DST  and  FUDS  datasets  of  the  eight  LiFeP04  cells  are  shown  in 
Figs.  19-21,  respectively.  The  first-order  RC  model  with  one-state 
hysteresis  was  again  found  to  be  the  preferred  choice. 

6.  Conclusions 

A  systematic  comparative  study  of  twelve  lumped  battery 
models  is  conducted.  The  MPSO  algorithm  is  applied  to  conduct 
parameter  optimizations.  The  model  comparison  is  performed  for 
two  types  of  Li-ion  cells,  to  assess  their  accuracies  in  both  training 
and  validation  datasets,  and  generalization  to  multiple  cells.  Cells  of 
the  same  chemistry  seem  to  have  very  similar  behavior,  especially 
when  they  are  new.  It  is  thus  possible  to  generalize  the  model  estab¬ 
lished  for  a  single  cell  to  predict  the  behavior  of  several  adjacent 
cells.  The  ability  for  a  model  to  be  generalized  is  thus  important. 
Since  complex  models  typically  are  more  expensive  and  more  sus¬ 
ceptible  to  uncertainties,  models  that  are  accurate  enough  and  yet 
simple  are  preferred.  Comparison  results  indicate  that  the  first- 
order  RC  model  is  preferred  for  LiNMC  cells,  while  the  first-order 
RC  model  with  one-state  hysteresis  seems  to  be  the  best  choice  for 
LiFeP04  cells.  The  developed  cell  voltage  models  can  be  used,  e.g., 
in  SOC  estimation  in  battery  management  systems. 
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